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ABSTRACT 

Infrared photometry and spectroscopy covering a time span of a quarter cen- 
tury are presented for HD 31648 (MWC 480) and HD 163296 (MWC 275). Both 
are isolated Herbig Ae stars that exhibit signs of active accretion, including driv- 
ing bipolar flows with embedded Herbig-Haro (HH) objects. HD 163296 was 
found to be relatively quiescent photometrically in its inner disk region, with the 
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exception of a major increase in emitted flux in a broad wavelength region cen- 
tered near 3 /jm in 2002. In contrast, HD 31648 has exhibited sporadic changes 
in the entire 3-13 fim region throughout this span of time. In both stars the 
changes in the 1-5 j^m flux indicate structural changes in the region of the disk 
near the dust sublimation zone, possibly causing its distance from the star to 
vary with time. Repeated thermal cycling through this region will result in the 
preferential survival of large grains, and an increase in the degree of crystallinity. 
The variability observed in these objects has important consequences for the in- 
terpretation of other types of observations. For example, source variability will 
compromise models based on interferometry measurements unless the interfer- 
ometry observations are accompanied by nearly-simultaneous photometric data. 

Subject headings: circumstellar matter - infrared:stars - planetary systems: pro- 
toplanetary disks - stars:formation - stars: pre-main-sequence - starsdndividual 
(HD 31648, HD 163296) 


1. Introduction 

HD 31648 (MWC 480) and HD 163296 (MWC 275) are “isolated” Herbig Ae/Be stars 
(HAeBes) that bridge the gap between classical pre-main sequence (PMS) and main sequence 
(MS) stars in the sense that, although they are Ae stars that lie above the MS and exhibit 
emission lines and infrared excesses, they are not immediately associated with dark clouds 
or reflection nebulae. 

With a distance of 1221^ pc and an age of 4^.5 Myr (van den Ancker et al. 1998), HD 
163296 appears to be an outlying member of the Upper Cen-Sco OB association. Its spectral 
energy distribution (SED) at wavelengths shortward of 15 /am is very similar to that of the 
classical HAeBe star AB Aur (Sitko 1981) with a strong continuum infrared excess and 10- /mi 
silicate emission band. HD 163296 is often considered a prototype of the “isolated” HAeBe 
star class. At longer wavelengths, the SEDs differ: HD 163296 lacks the extended “plateau” 
exhibited by AB Aur between A ~ 20 - 60 /on (Meeus et al. 2001). Recent coronagraphic 
observations show the presence of an underlying disk of dust, and outflowing gaseous jets 
in Lya (Devine et al. 2000) that include a string of Herbig-Haro (HH) objects (Grady et al. 
2000). At visible wavelengths, the dusty disk can be traced to 450 AU (3.7 arcsec) (Grady 
et al. 2000). Its inclination has has been reported to be between 45° and 65° (Mannings & 
Sargent 1997; Grady et al. 2000; Wassell et al. 2006; Isella et al. 2007), and we have adopted 
a value of 50° for our modeling. It also exhibits both far-UV emission lines (Deleuil et al. 
2005) and X-ray emission that have been attributed to optically-thin shock-heated gas that 
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is accreting onto the stellar surface, most likely along magnetic field lines (Swartz et al. 
2005). 

HD 31648 has a distance based on Hipparcos measurements of pc and age of 

2.5i('o Myr (van den Ancker et al. 1998), which would place it at the edge of the Tau-Aur T 
association. Using this same distance, Blondel & Djin A Djie (2006) derive an age of 9 Myr. 
However, Simon et al. (2001) suggest a distance of 170 pc and an age of ~7 Myr. It also 
exhibits an outflow of HH objects (Stecklum et al. 2007) and emission lines of Si IV (Sitko 
et al. 1981) consistent with hot accreting gas. Its dusty disk has recently been detected 
in scattered visible light (Grady et al. 2008), but it is faint and has a dropoff in surface 
brightness with radial distance that is steeper than most systems detected so far. The dust 
continuum is also detected 1.3 mm (Mannings & Sargent 1997). 

For both objects, the probable ages correspond roughly to the time when grain growth 
into larger bodies (i.e. cometesimals) has been operating, and in some models (Boss 1997) 
planetary embryos and perhaps true planets have already formed. 

Analysis of the infrared SEDs of Herbig Ae stars obtained with the Infrared Space Ob- 
servatory (ISO) show that most of these fall into one of two broad groups - those whose 
infrared emission that can be fitted (approximately) with simple power-laws, and those that 
require the addition of another roughly “blackbody” component (Meeus et al. 2001; Acke &; 
van den Ancker 2004). These form the Group II and Group I sources in the nomenclature 
of Meeus et al. (2001). Dullemond & Dominik (2004) suggest that these represent an evo- 
lutionary scenario, where Group II objects began originally as Group I objects with flared 
disks irradiated by the star, but as the disks age, grain growth decreases the optical depth 
of the disk, which then flattens and is eventually shadowed by a “puffed- up” inner rim of 
the disk. This inflated inner disk wall results when a paucity of accreting gas inside the dust 
sublimation zone (hereafter DSZ) leads to direct illumination of the dust by the star, causing 
the location of the transition zone to migrate outward and to expand vertically, compared 
to one that is shadowed by the accreting gas [see Fig. 2 of Millan-Gabet et al. (2007) for 
example] . 

Eventually, inner disk clearing, which can occur beyond the dust sublimation radius, will 
give rise to a third class of SEDs: those where the inner region is cleared due to the cessation 
of accretion, as well as the possible clearing by objects of planetary mass. Calvet et al. (2002) 
and Bergin et al. (2004) have recently provided strong evidence for the clearing of the inner 
disks of a number of PMS stars. In the disks surrounding TW Hya, GM Aur, DM Tau, 
and LkCal5, the spectral energy distributions require substantial clearing of material for 
distances closer to the star than a few AU. This process gives rise to the recently-recognized 
“transition disk” objects. Grady et al. (2007) find similar evidence for the star HD 169142, 
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another such transition object with a substantial dust wall near 30 AU, and an inner disk 
with a tremendous paucity of dust. In fact, many objects in the original Meeus Group I are 
not young objects with flared disks, but actually stars with transition disks with thick cool 
disks with large inner radii. 

The SEDs of HD 31648 and HD 163296 are consistent with them being members of the 
Meeus Group II. HD 163296 had the faintest disk of all of the Ae stars in the sample of 
Grady et al. (2005), with a middle zone that is darker than the outermost traceable reaches 
of the disk, possibly signaling its partial shadowing by the inner disk. The faintness and 
radial dependence of the surface brightness in the disk of HD 31648 (Grady et al. 2008) are 
also consistent with significant dust settling and shadowing. 

It is surprising, then, that both HD 31648 and HD 163296 exhibit many indicators 
of the active accretion generally associated with flared unshadowed Group I objects: hot 
non-stellar gas, jets, and HH objects. In fact, of the 4 HAeBes with known HH knots (HD 
163296, HD 31648, HD 104237, and AB Aur), the first three are Group II sources, in seeming 
contradiction to the general concept hat these Group II objects result from very low accretion 
rates. Apparently, the accreting gas can in some instances does not shadow the inner dust 
wall sufficiently, so that gas accretion and shadowing are not mutually exclusive. This might 
occur if the inner gas disk has a small scale height compared to the DSZ, or if it becomes 
optically thin in the radial direction. The degree to which the inner gaseous accretion disk 
can actually shadow the dust further out is not known precisely, but Muzerolle et al. (2004) 
suggest that the gaseous disks becomes optically thick radially once the accretion rate drops 
below a few times 10 -8 M 0 yr -1 . Its vertical scale height is sufficient to shadow the dust 
disk midplane, but not the entire inner wall of the dust disk. 

Highly collimated outflows require an efficient collimating mechanism. In the T Tauri 
stars this mechanism is generally thought to result from magnetic fields in the star, possibly 
connected to the disk (Shu et al. 1994). Hubrig et al. (2006) reported the detection of 
circular polarization in the circumstellar Ca II lines, but the implied derived field strength 
in this outflowing material is uncertain, as the polarization detection was only at a 2-cr level. 
No photospheric magnetic field was detected. They reported a 4-<r detection of a magnetic 
field in HD 31648, also mostly of circumstellar origin. However, Wade et al. (2007), with 
comparable observational uncertainties, found no detectable field. Only further study will 
reveal if these stars have significant, but variable, fields. 

The spacing and proper motions of the HH knots in HD 163296 (Grady et al. 2000; 
Swartz et al. 2005) requires changes in ejection activity on time scales of about 5-6 years. 
If these are controlled by the inner disks (inner gas disk or the gas+dust disk further out), 
then changes in disk structure may fuel the outbursts. On the other hand, if the outbursts 
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are driven by the physics of a magnetic field (such as reconnection events) rather than the 
fuel supply, then outburst activity might affect the inner disk structure through increased 
far-UV and X-ray activity and changes in disk winds generated by the activity. Either way, 
one might expect changes in the innermost regions of the disk in such objects on time scales 
that are similar to that of the generation of new HH objects. 

As part of a long-term program to investigate the mid-IR (~3-13pm) brightness vari- 
ations of HAeBe and related stars, we have observed HD 31648 and HD 163296 at these 
wavelengths on multiple epochs. Here we describe the nature of the observed variations, and 
how these relate to disk structure and accretion activity. The complete sample of objects 
observed as part of this program will be discussed in a later paper. 


2. Observations 

Observations were carried out over a span of a quarter century using a variety of pho- 
tometric and spectroscopic instruments. Table 1 lists the dates of the observations and the 
instruments used in each case. 

The bulk of the data presented here that cover the 3-13 pm region were obtained between 
1994 and 2006 using the Aerospace Corporation’s Broad-band Array Spectrograph System 
(BASS) on NASA’s Infrared Telescope Facility (IRTF). The 2004 BASS data were obtained 
at the Mount Lemmon Observing Facility (MLOF). This instrument uses a cold beamsplitter 
to separate the light into two separate wavelength regimes. The short-wavelength beam 
includes light from 2.9-6 pm, while the long- wavelength beam covers 6-13.5 pm. Each beam 
is dispersed onto a 58-element Blocked Impurity Band (BIB) linear array, thus allowing for 
simultaneous coverage of the spectrum from 2.9-13.5 pm. The spectral resolution R = A/AA 
is wavelength-dependent, ranging from about 30 to 125 over each of the two wavelength 
regions. The simultaneous coverage of the entire 3-13 /xm region that BASS provides has 
played a critical role in the investigations of the clearing of the inner disks in systems like 
those we are discussing here (Calvet et al. 2002; Bergin et al. 2004). 

At the IRTF, the circular entrance aperture of the instrument subtends 3.4 arcsec on 
the sky (8.5 arcsec at the MLOF) . At the distance of HD 163296 the BASS/IRTF 1.7- 
arcsec radius corresponds to approximately 210 AU, while the disk can be traced to 450 
AU in coronagraphic images (Grady et al. 2000). However, the disk has a half-width at 
half-maximum brightness of 0.007 arcsec at 12.5 pm (Leinert et al. 2004), and a half- width 
half-maximum of less than 0.5 arcsec at 1.3 mm wavelength (Mannings & Sargent 1997). 
In this paper we will limit our discussions to wavelengths of 20 pm and shorter, since the 
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entrance aperture includes essentially all the radiation emitted at wavelengths shortward 
of 20 fJL m. A wavelength of 20 qm corresponds (in units of A F\ Wm 2 ) to a characteristic 
blackbody temperature of T ~ 180 A", the equilibrium temperature at 12 AU (0.1 arcsec for 
a distance of 122 pc) from a star with L ~ 2 6L Q (van den Ancker et al. 1998). 

Both the 1.3 mm interferometry (Mannings k. Sargent 1997) and coronagraphic scattered 
light imaging (Grady et al. 2000; Stecklum et al. 2007) of HD 31648 indicate that it has an 
even smaller disk than HD 163296. Because all of the other entrance apertures used in this 
study were larger than this value, the flux at the wavelengths sampled for variability (i.e. , 
less than 20 /jm), will be contained within all of the entrance apertures used. 

Prior to 1996, the data on both stars were obtained using primarily medium-bandpass 
filters with spectral resolutions of R ~ 6, and narrower-bandpass filters with R ~ 50 near 
the 10-qm silicate band. When searching for possible variability in the infrared emission, 
it is critically important to insure that all the observations are based on the same relative 
and absolute flux systems. Unfortunately, throughout the history of infrared astronomy, a 
wide variety of filters, standard star magnitude systems, and absolute flux calibrations have 
been in use. For this study, we have applied a single relative brightness and absolute flux 
calibration system, described more fully in Appendix A. 

High (R~1000) spectral resolution measurements at shorter wavelengths were carried 
out in 2002 and 2005 using The Aerospace Corporation’s Visible and Near-Infrared Imag- 
ing Spectrograph (VNIRIS). The original infrared portion of the spectrograph (NIRIS) is 
described by Rudy et al. (1999). It employs two independent channels separated by a beam- 
splitter covering the wavelength ranges 0.8-1. 4 /u,m and 1.4-2. 5 /im, respectively. Each channel 
employs its own collimator, camera, grating, and Hawaii 1 HgCdTe detector array, but views 
the slit through a common field lens. This was the configuration for the 2002 observations 
of HD 163296. By 2005, a second beamsplitter and a third channel, covering the 0.4-0. 9 
li m, had been added to the system, now called VNIRIS (a term which we will use for both 
versions throughout the remainder of this paper), and used to observe the same star that 
year. It resides outside of the cryogenic vessel but also views the same slit through the same 
field lens as the infrared channels. The detector in this case is a red-sensitive, deep-depletion 
CCD. 

Dispersion is nearly uniform over the wavelengths covered by each channel of the VNIRIS 
spectrograph, implying that resolution the R = A/AA increases linearly with wavelength. 
For objects that underfill the slit width, resolution is determined by seeing and guiding. 
We measured the resolution directly using the unresolved emission lines of HD 163296. 
For the July 2002 observations it ranged from approximately 550 to 950 for the 0.8-1. 4 fim 
wavelength range spanned by the short wavelength channel; for the long wavelength channel, 
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the resolution varied from about 500 at 1.4 pm to 900 at 2.5 p m. For the 2005 observations 
the resolution for the infrared channels was very similar to that in 2002. For the optical, 
which did not exist for the July 2002 observations, the resolution ranged from 600 to 1000 
between 0.5 pm and 0.9 pm. Data shortward of 0.5 pm were acquired in second order; the 
resolution varied from 900 to 1150. 

The July 2005 observations used the solar type star HD 162255 for a standard. For 
a model of its spectrum, we used that of the sun and modified the shapes and strengths 
of the stronger absorption lines (e.g., Ha, the infrared Ca II triplet) to match the actual 
observations of the standard. The same process was applied to the earlier type standard 
used for the July 2002 observations, but for it the main focus was the H I lines. 


3. Results 

The SEDs of HD 163296 and HD 31648, from ultraviolet to millimeter wavelengths, are 
illustrated in Figures 1 and 2. For both objects, the flux in units of AF a Wm ~ 2 is lower 
at 60 pm than at 8 pm, consistent with a Group II classification. Neither object exhibits 
significant variation in the photometric V band at 0.55 pm: historically these have been 
restricted to being less than 20% for HD 31648 and 10% for HD 163296 (de Winter et al. 
2001). In the near- infrared (JHK bands), the variations are larger for HD 163296 than in 
the V band, being over 30% at J(1.2 pm), 20% at H(1.6 pm), and 15% at K(2.2 pm) (de 
Winter et al. 2001). HD 31648 has not been extensively monitored, and little information 
is available from the literature on its near-IR variability, but based on our limited data set, 
HD 31648 seems to exhibit near-IR variability at least as large as that of HD 163296. 


3.1. HD 163296 

For HD 163296, Figure 3 shows data obtained on 14 dates, where we have supplemented 
our observations with select data from the literature. Many sets of these observations were 
obtained within a few weeks of one another, and together could be considered single “epochs” , 
providing coverage on a total of 10 epochs. Observations by our team began in 1979 with the 
Kitt Peak bolometer, and after a hiatus of 17 years, resumed with BASS in 1996. Beginning 
in 1996, we obtained BASS spectra on 6 separate epochs. On two of these, 2002 and 2005, 
we also observed the star with VNIRIS within three weeks of the BASS observations. In 
addition, the Carlos Sanchez telescope of the Teide Observatory on Tenerife obtained JHK 
data (Eiroa et al. 2001) on this star within 4 days of the BASS observations in 1998. The 
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Infrared Space Observatory (ISO) Short Wavelength Spectrograph (SWS) data in Figure 
3 were obtained within 4 days of the 1996 BASS spectrum as well, and are in very good 
agreement at all wavelengths in common. We have also included three epochs of JHKLM 
photometry from de Winter et al. (2001) between 1980 and 1986. 

Most of the data are consistent with no variability in excess of 10%. In 2002, however, 
there was a very significant rise in the flux level between 1 and 5 /am, with perhaps only 
a slight increase in flux in the 10 pm silicate feature. We were fortunate to have VNIRIS 
and BASS observations in both low-flux (2005) and high-flux states (2002). The combined 
VNIRIS-BASS data for the two epochs are shown in Figure 4. No shifts in flux level has 
been applied to these observations in order to make them “match” at any wavelength. The 
joint VNIRIS-BASS spectra clearly show that the “outburst” in 2002 was seen in both 
instruments at both telescopes, and is unlikely to be the result of instrumental effects and 
data processing. 

The fairly regular spacing of the HH objects associated with the outflow jet in HD 163296 
requires quasi-periodic ejection events roughly every 5-6 years, de Winter et al. (2001) report 
JHKLM photometry during three epochs in the gap on our own coverage between 1979 and 
1996 (labeled as “ESO” in Figure 3). Figure 5 illustrates the changing brightness of the 
1-5 pm emission as measured at 3.77 pm using the ESO and BASS data (the 1979 KPNO 
observations were at 3.5 pm and of lower quality, and have been omitted here). Two of 
the ESO points (1980 and 1983) are clearly consistent with the low flux state seen in most 
epochs, while their 1986 observations seem to have been obtained in another earlier outburst 
state. The data seems consistent with a variability time scale of about 16 years, although 
shorter time scales cannot be ruled out due to the paucity of the time coverage. There is 
also insufficient temporal coverage to decide if the changes are periodic. Hence, we cannot 
yet tell for sure if the variations observed are related to the quasi-periodic ejection of the 
HH objects. 


3.2. HD 31648 

For HD 31648, 10 epochs of data are shown in Figure 6. Our temporal coverage for 
HD 31648 is less extensive than for HD 163296. We include JHK photometry from the 
2MASS survey and data obtained with the ISO SWS. Of the three sets of BASS data that 
simultaneously cover the 3-13 pm region, none shows evidence of a truly strong, well-defined 
1-5 pm hump (except perhaps in its highest flux state) that is the hallmark of objects whose 
SEDs in that spectral region are interpreted as being due to the irradiated puffed-up inner 
disk wall, and which is usually used to invoke the self-shadowing picture of the Meeus Group 
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II sources. As will be discussed below, the weakness of this feature may also be inconsistent 
with the interpretation of visible-wavelength interferometry measurements, which generally 
predict a substantially stronger feature (Monnier et al. 2006; Isella et al. 2006). 

Unlike HD 163296, there is yet no indication of a “preferred” lower flux state. The 
observed variability in HD 31648 is also of greater relative amplitude than that observed so 
far in HD 163296, and is more prominent in the 10 pm. silicate band, as shown in Figure 
4. Figure 7 shows the BASS-derived fluxes at 3.77 pm and 10.5 /iin for both stars. For HD 
31648, the flux in the 3 yum and 10 yum regions vary together, whereas the only significant 
change in the 10 pm emission of HD 163296 occurred in the 2002 outburst. Nevertheless, the 
similarity in the “difference” spectra indicate that we are seeing the same basic phenomenon 
on both objects. 


3.3. Monte Carlo Radiative Transfer Models of the Disks 

In order to better understand the nature of the dusty disks surrounding HD 31648 
and HD 163296, we have used of the Monte Carlo radiative transfer code of Whitney et 
al. (2003a, b, 2004) to model the SEDs of these systems. The code is a simulation using 
many photons propagating through, and interacting with, a circumstellar environment to 
produce the SED. The circumstellar dust model is separated into three physical zones: a 
disk, an envelope, and an outflow region. Each region is assigned user-defined geometrical 
dimensions and grain properties. For the disk, the inner and outer radii, inner edge disk 
scale height, mass, and the amount of flaring can be set. When hydrostatic equilibrium and 
mass conservation in the accretion flow are enforced, the radial mass density and scale height 
are coupled. In the case of an irradiated a-type disk, the surface density E oc R~ l (D’Alessio 
et al. 1998) where R is the radial distance. For a radial mass density p oc R.~ A and density 
scale height H oc R B , mass conservation in the accretion flow requires A — B = 1 (Robitaille 
et al. 2006), a restriction that was imposed in these models (see Table 2). Changing B 
alters the degree of flaring, while maintaining the restrictions of hydrostatic equilibrium and 
conservation of mass accretion at any given radial distance. B = 1 describes an unflared disk 
(i.e. conical opening), while B > 1 is for a standard “concave upward” flared disk. B < 1 
can be used for an “anti-flared” convex upward disk. The last of these may be appropriate 
for disks with grain settling. Envelope properties include inner and outer radii and the rate 
of flow of material onto the star. Such an envelope can also be used to simulate a disk wind 
or mini-halo. The outflow region is useful for including the effects of jets, as is observed 
in many young disk systems. This dust model is then overlaid on a radial/azimuthal grid, 
and the material density of each region is calculated for every grid cell. The main disk 
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itself is a blend of two regions, the disk midplane and an upper disk region analogous to a 
disk “atmosphere” discussed by many investigators. Transition between grain characteristics 
going from midplane to the disk surface can be implemented when the mass density crosses 
a specified threshold, which in our model of HD 163296 was H 2 density = 10 10 cm -3 , and 
larger sized dust grains were assigned to the disk midplane than the atmosphere to simulate 
the effects of dust settling. For HD 31648, the grain size was uniform throughout the disk. 

All of the dust property files used here consist of nearly half and half mixture of amor- 
phous carbon and silicates. Three dust files (of the four that come with the standard release 
of the radiative transfer code) with differing grain size distributions are used in the model- 
ing, although it is not necessary (or desirable) to use them all in every object studied. The 
first is the interstellar medium (ISM) grain model of Kim et al. (1994) (hereafter KMH), 
which utilizes a power law size distribution with an exponential cutoff. The second one, 
built on the KMH model, is the one used by Cotera et al. (2001) for modeling the scattered 
light in the upper disk atmosphere of HH 30. This one has a flatter size distribution, and a 
mean grain size about twice that of KMH - the size where the distribution has a turnover 
is near 0.55 pm, as opposed to ~ 0.2 pm for KMH. Finally, a “large” grain model of Wood 
et al. (2002) for the infrared SED of HH 30 was used, which is similar to the other two, but 
with the turnover near 50 pm. In each region, the grain temperature is determined by the 
locally-generated accretion, irradiation from the star, and re-radiation from within the disk. 
For simplicity, a single temperature is used in each cell, based on the mean grain properties 
(Bjorkman & Wood 2001). In the fitting process, the SED in the millimeter-submillimeter 
will govern the use of smaller versus larger grains. Figure 8 shows the structure of the inner 
disk of the model for HD 163296, with the grain size regimes illustrated. 

The code also allows for the presence of accretional heating of the star. This includes 
both a magnetically channeled accretion column, and collisional heating of the stellar pho- 
tosphere. These are modeled using the ideas of Calvet & Gullbring (1998). The accretion 
column produces X-rays, which for simplicity is assumed to radiate uniformly at 100-500 
A over an area a fraction of that of the star, while the heated atmosphere is emitted as a 
blackbody uniformly over the star. 1 

Many of the model input parameters use for these two sources are restricted using 
imaging information. The outer radius of the disk was determined using coronagraphic 
imaging obtained using the Hubble Space Telescope (Grady et al. 2000; Stecklum et al. 
2007). In the case of HD 163296, the scattered light images reveal a radial surface brightness 
distribution that is somewhat radially-dependent and azimuthally-dependent, but close to an 


1 described in the code updates: http://gemelli.colorado.edu/~bwhitney/codes/codes.html 
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R~ 3 power law. This requires a disk scale heigh parameter B ~ 1.0. For our model, B = 0.99, 
and predicts the correct surface brightness distribution at radial distances sampled by the 
imaging. The disk of HD 31648 was undetected in the survey of Grady et al. (2000), but later 
detected by Stecklum et al. (2007) and found to have a steeper radial dependence consistent 
with greater dust settling. For HD 31648, our model uses B = 0.6, an “anti-flared” disk. In 
the case of HD 31648 the predicted surface brightness is not as steep as the observed one, 
and may be affected by the presence of grain scattering phase functions that are different 
than those of the model grains. Because grain growth and coagulation are likely to lead to 
porous non-spherical grains, investigating the properties of those grains in the context of the 
radiative transfer in disks is warranted. 

Similarly, the disk inclination for HD 163296 was set to be 50°, within the range (45° 
and 65°) reported in the literature (Mannings &; Sargent 1997; Grady et al. 2000; Wassell et 
al. 2006; Isella et al. 2007). For HD 31648, the 30° inclination is close to the 38° of Simon 
et al. (2001) and 31° of Pietu et al. (2006) (both with quoted uncertainties of ±1°). 

In the models for HD 31648 and HD 163296, the model outflow regions are empty of 
material (we are not attempting to model the jet material) , but a small envelope has been 
included, with an outer radius that is much smaller than that of the disk. In these models, 
because these two stars are “isolated” from their parent star-forming clouds, we do not 
expect a true “infalling envelope” of material that this component of the model is usually 
used to describe. Here, we are simply using it to simulate the presence of material outside 
the disk in the inner few AU, where the variability activity is occurring. The inner disk 
edge is directly illuminated by stellar light and therefore creates excess flux in the 1-5 /mi 
region, which is seen in the SEDs of both stars (albeit to a much smaller degree in HD 31648 
than in HD 163296). The small envelope provides some of the flux in the 8-12 gm region. 
The emission from this envelope also reheats the portions of the disk farther from the star, 
creating an excess of disk flux over that from stellar emission alone in the mid and far IR 
regions. This component was necessary for obtaining adequate fits in the near-mid IR, and 
its existence, at least in the case of HD 163296, is also indicated by the extended “halo” 
component detected by the interferometry observations of Monnier et al. (2006). However, 
in all of the models, the disk surface dominates the flux in the 3-13 pm region. In the case 
of HD 163296, for example, 90% of the flux is from the disk, and 10% from the envelope. For 
the disk flux, less than 1% is from accretional heating. The rest is from direct irradiation 
from the star, supplemented by additional warming from the envelope. Figure 9 illustrates 
the temperature and density structure for the HD 163296 model. 

It should be noted that in these models the disk masses are higher than those usually 
derived from millimeter data. For example, Mannings et al. (1997) provide two estimates 
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of the disk mass for HD 31648. The first assumes an optically thin source at 1.3 mm, a 
radially-averaged temperature of 40 K, and a grain opacity k v of 0.1(zv (GHz)/ 1,200)^ cm 2 
g -1 . This yields a mass of gas and dust of M ~ F u d 2 / k u B w (T) = O.O24M 0 for a gas/dust 
ratio of 100 and 0=1. A second estimate using power law radial temperature and density 
gradients results in a disk mass about twice as large. These estimates depend on the assumed 
dust opacity law k v (among other things). In the case of HD 31648 the grain file used for 
our models has a 1.3 mm opacity only 3% that used by Mannings et al., which results in 
a substantially more massive disk. Use of the “large” grain population for the disk mid 
plane (as in HD 163296) would lower the disk mass slightly, but not substantially, since the 
emission is dominated by the warm disk surface in these models. Extending a large grain 
model to the surface would result in a 10 //m band much weaker than that observed, unless 
the “envelope” were made more massive. 

In this paper we investigate two possible sources of NIR variability in these systems. 
These are similar to those discussed by Vinkovic & Jurkic (2007) for low luminosity systems, 
and are illustrated in their Figure 12. In their model for low-luminosity systems, they 
postulate both (A) a puffed-up inner disk and (B) the ablating of a portion of the dust in 
the inner disk region into a halo above the disk. We will show that changes in one or both 
may be the source of the variations seen. 


3.3.1. Changing Dust Sublimation Zone (DSZ) 

In the case of HD 31648, we are dealing with an SED that is well into the Meeus Group 
II category, and generally consistent with the faintness of the outer disk in scattered light. To 
meet these restrictions we used a disk that was essentially anti-flared, its surface curvature 
defined by a surface scale height h which varies a radial distance r as H oc R 0 6 (slim disk) 
and extends from 0.2 [smaller than that measured interferometrically by Isella et al. (2006)] 
to 250 AU. For acceptable fits to the size of the 10-yum silicate band and overall SED, the 
entire disk utilized the grain model of Cotera et al. (2001), as using the larger grains of Wood 
et al. (2002) led to a poorer fit at the longest wavelengths. The envelope above and below 
the disk extends from 0.2 to 10 AU and uses the ISM-like grains of KMH. The input star 
temperature was 8250 K. Accretional heating of the star was included. This heated stellar 
photosphere produces much of the flux in the 0. 1-0.2 /am region, and emission of this nature 
is required for matching the UV flux in this star (Blondel & Djin A Djie 2006). For HD 
31648, the “coverage” of the X-ray accretion column was 4%. 

Because the details of the fit of the inner region of the disk will depend on the data set 
used, in Figure 10 we show the model for the 1996 BASS observations only. In Figure 11 we 
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show both the 1996 and 2004 BASS observations, along with their respective models. The 
model for the epoch 2004 data is essentially the same as that for 1996, but increasing the 
scale height of the disk at its inner edge from 0.0076 AU to 0.013 AU, and a slight increase 
in material in the envelope. These changes were introduce solely to be able to approximate 
the change in the SED. Changes in the inner structure of the disks in some PMS stars has 
been detected in recent interferometric observations (Tannirkulam, private communication), 
but simultaneous interferometry /SED observations will be needed for confirming the precise 
nature of these changes. It is currently unclear if there was a significant change in the 
location of the inner radius (through increased grain sublimation) , because the weakness of 
the NIR hump compared to the star’s photospheric emission makes constraining the radius 
of the inner disk difficult using SED information alone. 

In the current version of the radiative transfer code, the disk scale height is set by the 
requirement that the disk be in hydrostatic equilibrium throughout the disk, and setting 
the density scale height at the inner edge determines the scale height throughout the entire 
disk since it uses a single scale height power law exponent. In reality, the sound speed is 
insufficient for the entire disk to respond so quickly, and the system cannot be characterized 
by a single scale height power law. Models that include such dynamic processes are beyond 
the scope of this investigation. 

It is also unrealistic to attempt to fit the entire SED in complete detail in any case, as the 
lack of simultaneous observational coverage at all wavelengths would make such comparisons 
of limited value, except to derive disk characteristics in the most general way. In principle, an 
increase in scale height of the inner disk wall should alter the degree of shadowing of the outer 
disk, and hence its emission. Obtaining complete simultaneous 1-100 n m SED coverage will 
be difficult. Coordinated 1-13 /urn spectral coverage coupled with coronagraphic observations 
of the scattered light in the outer disk, however, may be a viable alternative to test for the 
presence of variable shadowing. 

For HD 163296, the disk is marginally anti-flared, its surface curvature defined by H oc 
R° ", and the disk extends from roughly 0.3 AU [compared to 0.45 AU derived from the 
interferometry of Monnier et al. (2006)] to 450 AU. The envelope extends from 0.9 to 10 
AU and consists of KMH grains. The disk begins with KMH grains near its surface, but 
transitions to the large grains of Wood et al. for the midplane [use of the smaller grains 
of Cotera et al. (2001) alters the emission only slightly at sub- millimeter wavelengths]. For 
the model fitting, we used the “quiescent” state BASS spectrum from 2005. The input star 
temperature was 8750 K, consistent with the spectral type of the star. The accretional 
heating near the star is less noticeable in this star, due to its earlier spectral type and 
subsequently greater intrinsic photospheric UV flux. This model is shown in Figure 12. The 
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2002 “outburst” can be fit with the same model, but with an increase in disk height from 
0.012 AU for the disk state in 2005 to 0.018 AU for the earlier 2002 state (Figure 13). In 
the case of HD 163296, the greater strength of the 1-5 /am emission makes the location of 
the inner dust sublimation zone somewhat easier to constrain than in HD 31648, and was 
different in the two epochs, being roughly 0.29 AU in the low flux stare and 0.35 in the 
higher state. The fact that these values are some what smaller than those of Monnier et al. 
(2006) probably indicate that part of the shortest-wavelength emission of the 3 /iin hump 
is contaminated by emission from hot gas interior to the DSZ. This gas has been recently 
resolved using interferometers (Tatulli et al. 2007; Tannirkulam et al. 2007). 


3.3.2. Detachment of Dusty Material in a Halo 

In the previous section we assumed that a change in the location of the DSZ produced 
the observed variation in the SED. However, such variations might also result from dust that 
is removed from the inner wall region and driven outward, as suggested by Vinkovic & Jurkic 
(2007) in their component (B). In Figure 4 we show the difference in the flux levels in HD 
163296 between 2002 and 2005, and model the expected thermal emission of grains capable of 
producing such a difference, using a simplistic “model” based on standard Mie calculations 
for the expected emissivity (Bohren & Huffman 1983). The net emissivity will depend 
on the size distribution of the grains and their mineralogical content. For our model we 
have used the optical constants of Henning et al. (1999) for graphite and pyroxene silicates, 
( Mg x Fei- x )Si0 3 with x = 0.8. Because the disk has probably entered the phase where grain 
growth and sticking may have proceeded onto the cometesimal-forming phase, we adopted 
a particle size distribution consistent with those measured in situ for Comet lP/Halley 
(McDonnell et al. 1987, 1991). The maximum and minimum size limits (which were assumed 
for simplicity to be the same for both materials), the silicate/graphite abundance ratio, and 
the (single) temperature of the grains were adjusted to reproduce as best as possible the 
general shape of the observed difference. The final grain size ranged from 0.03 /rm to 25 /im, 
with a silicate/graphite ratio by number of 0.5 and a grain temperature of 800 K. Our goal is 
not to create a definitive description of the “outburst” component of the SED, but simply to 
illustrate that the major features are consistent with the presence of an added population of 
small grains illuminated by the star. Such grains, when driven away from the star by either 
winds or radiation pressure, may account for the “halo” components suggested by Monnier 
et al. (2006) for some PMS disk systems, which in the specific case of HD 163296 was 5.0% 
± 2.5% at A=1.65 firn. Also shown in Figure 4 is a similar model for the difference in flux 
for HD 31648 between 1996 and 2004, based solely on the 3-13 /am BASS data. 
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For comparison with the model for HD 163296 described in the previous section, both 
models are shown close up in Figure 13, while in Figure 14 we show the entire SED for 2002, 
with the “outburst” grains providing the extra emission in 2002. 


4. Discussion 
4.1. HD 163296 

TV: ’ p,m spectral region that underwent the outburst in 2002 corresponds to the 

“hump” that is commonly seen in the SEDs of many HAeBe stars, and which is generally 
attributed to emission from the innermost “wall” of the dust disk (Natta et al. 2001), where 
the dust is subliming (the DSZ) and the incident stellar flux hits the dust layer orthogonally 
(Duliemond . Dominik 2004). If so, the brightening seen in 2002 would require either 
increased illumination from the star, or an increase in the surface area of the “wall” or some 
equivalent nearby structure. 

A comparison of the VNIRIS spectra in 2002 and 2005 indicate that the photospheric 
flux at 0.8 n m was constant, so that the outburst in the 1-5 /im region cannot be due to 
variable heating by the the star at visible wavelengths. However, because the accretion shock 
contributes to the heating of the disk, the possibility of heating by far-ultraviolet and X-ray 
photons produced in that region cannot be excluded. Otherwise, the variable emission must 
be due to changes in the disk itself, and in one instance we will discuss later (the T Tauri 
star DG Tau) only changes near the DSZ can account for the changes observed. 

Vinkovic & Jurkic (2007) discussed two ways that such variability might occur in low 
luminosity PMS disk systems. One is increasing the scale height of the inner edge of the 
inner dusty disk, which is puffed up because the low mass accretion rates result in the inner 
gaseous disk being optically thin in the radial direction. The other is by the removal of dusty 
material from the disk into a halo above the disk. Both are plausible, since the former is just a 
time- variable modification of the standard Duliemond & Dominik model for shadowed disks, 
which must occur if the shadowing by accreting gas changes, while the latter is consistent 
with the excess “halo” emission suggested by the, interferometry observations by Monnier et 
al. We have shown that both are capable of reproducing the general character of the near-IR 
variability, and both may be operating simultaneously. 

In both models, however, the presence of small silicate grains in the wall poses some 
problems for the standard puffed-up inner disk model. Recently, Vinkovic et al. (2006) have 
shown that a wall with a vertical inner edge can only explain the near-IR emission if the 
dust has perfectly gray opacity, which is contradicted by the presence of the silicate feature 
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at 10 /urn produced in this region, as it seems to be in these two objects. The characteristic 
temperature of the dust in the “outburst” in HD 163296 places the location of the grains 
at or very near that of the DSZ, and because the variations we see in HD 163296 occur 
simultaneously in both the hump and the silicate band, the simplest explanation is that the 
small grains are located at or near the wall, and that the opacity cannot be gray. As suggested 
by Millan-Gabet et al. (2007), this re-opens the debate concerning disks vs. envelopes as the 
source of the near-IR hump and other spectral features. Highly curved inner walls, which 
allow the incident radiation to arrive at a more glancing angle of incidence, may avoid these 
difficulties (Isella et al. 2006; Monnier et al. 2006), but until the models implied by the 
interferometry can reproduce the SED of this region more accurately, it is unclear whether 
invoking this geometry alone will solve this problem. 

In Figure 4 we apply the HD 163296 grain model to the difference between the 1996 
and 2004 flux levels for HD 31648. Here a higher temperature (T = 1600 K) was required to 
fit the 3-5 /mum SED. In this case, the grain model had a minimum grain size of 0.24 yum, 
as suggested by Isella, and a maximum size of 3 fim. It clearly underestimates the silicate 
band strength. If 0.2 yum were a “typical” grain size in this region, it would have to be the 
only grain size. Such a monodisperse grain population seems unlikely. 

Is the NIR variability related to accretion activity? Figure 15 shows the Ca II triplet 
lines near 0.85 /im, as well as the hydrogen Paschen (3 and Bracket 7 lines for HD 163296 
obtained with VNIRIS in 2002 and 2005. While these lines have usually been interpreted as 
a measure of accretion activity in PMS stars (Muzerolle et al. 1998, 2001), the net strength 
of the hydrogen bands did not change between 2002 and 2005, while that of Ca II were 
actually lower during the “hump outburst” of 2002, compared to the more quiescent state in 
2005. This would suggest that this measure of accretion activity can vary out of synch with 
increased disk wall emission. It has recently been suggested that much of the strength of 
the hydrogen emission lines may arise fairly far out in the gaseous disk (Tatulli et al. 2007; 
Tannirkulam et al. 2007), and be less coupled to changes in accretion activity very close to 
the star, where we would expect changes to be most rapid. 


4.2. HD 31648 

HD 31648 exhibits large variations in its IR flux over the 1-13 yum wavelength range, 
and its 1-5 yum hump implies a smaller DSZ scale height than in HD 163296. In the models 
of Dullemond & Dominik (2004), such a small hump, combined with an even weaker 60 yum 
flux, and with millimeter fluxes within 3.5 dex of the bump flux requires that most of the 
mass (99.9% or more) be in very large grains (although our model in Figure 10 did not use 
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the “large” grains of Wood et al. (2002) in the disk midplane). The faintness of the scattered 
light in the outer disk in the coronagraphic observations of Grady et al. (2008) is consistent 
the coagulation of small grains into large ones and a high degree of settling. In such a model 
the dust settling leads to a disk that is completely shadowed by the inner wall. 

The structure of the inner disk region of HD 31648 has been probed using interferometry 
at NIR wavelengths. Measurements obtained using the Palomar Testbed Interferometer in 
the K (2.2 jtim) photometric band (Isella et al. 2006), however, are only roughly consistent 
with the actual SED of this object. HD 31648 was also recently observed with the IOTA 
interferometer (Monnier et al. 2006) in the H photometric band (1.6 /im), and the observa- 
tions fit best using an inner scattering dust ring with rounded edge and diameter 0.50±0.04 
AU. However, this model also fails to reproduce the bulk of the SED photometric data on 
this object. When comparing the model SED with observed SEDs taken from the literature, 
Monnier et al. included no error bars in the literature SED data they plot, as the scatter 
in the SEDs (source variability) dominated over the individual measurement uncertainties 
of the SEDs themselves. This variability is clearly illustrated in our Figures 3 and 6, and 
reinforces the necessity of obtaining interferometry and photometric SED observations si- 
multaneously. However, as seen in these figures, it is not simply the object’s variability that 
poses potential difficulties for the models. At no time does HD 31648 exhibit a truly strong 
hump (spectral slope of zero near 2 pm in units of AF a ) except in its highest flux state 
(Figure 11). The success or failure of fitting the model SEDs derived from interferometric 
observations to photometrically-derived ones will depend on having contemporaneous pho- 
tometric data. Furthermore, a substantial amount of the emission at 2 pm emission seems to 
be coming from gas interior to the DSZ (Tannirkulam et al. 2007) in some of these systems, 
and this may apply to HD 31648. 

Other potential issues in model SEDs based on NIR interferometry of PMS disk systems 
have been pointed out by by Pontoppidan et al. (2007). One is that coverage of the u-v plane 
is not always optimal for the extraction of precise spatial scales. For example, using a number 
of different types of spatial models (disks, gaussians, rings, etc.), Eisner et al. (2004) found 
a range of a factor of two in acceptable angular scales for the “size” of the inner disk in HD 
31648. Isella et al. (2006) found a range of acceptable values for inner rim, 0.53-0.63 AU, for 
the same disk. 

Another issue is that of “grain survival”. Because smaller grains tend to be hotter 
than larger grains, there is likely a change in grain size with distance from the star, as 
larger grains can be closer without undergoing sublimation. In many of the models based 
on interferometry, only a single grain size is used. Isella et al. used a grain radius of 0.2-0. 3 
pm for their grain size, which they describe as “big” (i.e. the large end of ISM grain size 
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distribution). Pontoppidan use 0.3 pm grains for their small grains and 3.0 pm for their 
larger ones. They found that the distance to the inner rim of the disk in the Herbig star VV 
Ser derived by Isella et al. based on NIR interferometry was about half that derived from 
their own SED fitting. Their best-fitting VV Ser model for the SED and interferometry taken 
together has their big grains at 0.27 AU and small grains at 0.80 AU. A similar ambiguity 
may apply to HD 31648. In fact, the geometry used by Pontoppidan et al. shares many of 
the features of the models used here. 

Objects like HD 31648 present challenges to our understanding of disk structure and 
evolution. The weakness of the scattered light in coronagraphic imaging experiments and the 
overall shape of the SED would normally indicate significant dust coagulation and settling 
has occurred, features that are characteristic only of the later post-accretional evolutionary 
stages of disk evolution. Yet HD 31648 is still actively accreting and producing polar jets 
with HH objects. Any shadowing must also be accomplished with an inner dust sublimation 
wall of very small vertical scale, since the 1-5 pm hump in the SED is so weak. In our 
models it is less than 0.02 AU high. Taken together, these facts suggest that the disk in 
HD 31648 may be more like a “classical” accretion disk with an inner gas disk merging 
smoothly with the gas+dust zone, with the former providing enough shadowing to prevent 
the dust sublimation zone from becoming too puffed-up. But such actively accreting disks 
generally occur in young systems with highly flared disks, which is not the case here. A 
disk with the sort of extreme grain growth and settling that produces an SED like that HD 
31648 should no longer be actively accreting gas. Exactly how low the accretion rates can 
be and still provide adequate shadowing is addressed by Muzerolle et al. (2004), where the 
gas disk becomes radially optically thin in the mid-pane for accretion rates between 10~ 8 
and 1O~ 9 M 0 yr -1 . In our model the accretion rate is 4xl0~ 9 M 0 yr -1 (Stecklum et al. 2007), 
right where this transition may occur. Even in the model of Muzerolle et al., though, only 
a small fraction of the dust wall is actually shadowed. This star is likely evolving into a 
transition disk object. 


4.3. Disk Instabilities at the Dust Sublimation Zone 

As yet there is no clear explanation for the source of the changes in and around the 
DSZ that give rise to the variations on the emission observed in these systems. A number 
of possible sources of such changes are internal instabilities will be discussed in turn. 
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4-3.1. “Standard” a-Disk Instabilities 

Instabilities in long-lived disks can occur in a quasi-periodic manner if there is more 
than one “stable” equilibrium configuration that lead to a limit cycle behavior, as occurs in 
some cataclysmic variables. The most likely location for instabilities to occur is where there 
are significant changes in the radiative opacity, equation of state, or other important source 
parameters within the system. One such location is the transition from an optically thick 
region dominated by the dust opacity and the dust-free inner gaseous disk (which in principle 
might be optically thick or thin): the inner boundary of the DSZ. While it is beyond the scope 
of this paper to determine the precise physical mechanism for the behavior of these disks, 
we can place some useful constraints on the nature of these changes by comparing them to 
what is expected from different instability time scales. Because the variability is primarily in 
the region responsible for the NIR hump, the boundary between where the radiative opacity 
being controlled by dust to where it is being dominated by gas is implicated. 

Can instability-driven changes in the disk structure in the vicinity of the inner “wall” of 
the DSZ lead to the sort of changes observed? Because the disk surrounding HD 163296 is 
still actively accreting, it can be heated by both the radiation from the star and from internal 
viscosity. At the present time, the stability of disks heated both internally and externally is 
not well understood, but we can get an approximate idea of variability time scales, some of 
which do not depend explicitly on the degree of viscous energy dissipation. 

The general physics of such disks are discussed by Shakura & Sunyaev (1973), while 
Pringle (1981) outlines some of the instabilities and their timescales in these systems. Typ- 
ical timescales encountered in accretion disks are the dynamic timescale the hydrostatic 
timescale t z , the thermal timescate t t h, and the viscous timescale t u . The dynamic timescale 
is the time required for a parcel of mass in the disk (such as a hot spot) to move a significant 
fraction of its orbit, and is given by 

^ = 5 = n ‘ 1 

where R is the radius outward in the disk, v $ is the orbital velocity, and Q is the angular 
velocity (note that the orbital period is 27rt^). For R in AU, t# in years, and the stellar mass 
M in M©, Keplerian orbits will have 

+ - l R 3/2 

L< t > 2 tt M 1 / 2 ' 

Depending on the disk inclination and thickness, any hot spot might be more easily visible 
at some time than at others. 

The timescale required for a region of scale height H and sound speed c s to establish 
hydrostatic equilibrium is 



For changes in the 1-5 fi m feature, the relevant distance is the dust sublimation zone, which 
for blackbody-like grains will be located approximately at 

«~(¥) 2 ( 4) 1/2 AU - 

For a star with L ~ 2 6L 0 and a grain evaporation temperature of T ~ 1500 If , the inner 
dust wall would be at R ~ 0.3 AU. In this case, 

tfj, ~ t z ~ few days 

for a star whose mass is ~ 2.5 M 0 . 

By contrast, the thermal time scale is the time required for a disk to dissipate a sizable 
fraction of its internal heat, and is given by 

tth ~ ~t<!> 

where a is the standard dimensionless viscosity parameter (Shakura k. Sunyaev 1973). Hart- 
mann (2005) has suggested that, for a variety of reasons, a ~ 0.001 is a good reasonable 
global average for the accretion history in disks surrounding T Tauri stars, with 0.01 being 
an upper limit. If true for the systems being discussed here, a ~ 0.01-0.001 would give 

t t h ~ 10-100 years, 

Finally, the viscous time scale, over which the disk surface density can change, is 
U ~ (§ ) 2 tth- 

For a thin disk tv must be considerably longer (e.g., thousands of years). 

It is striking that the IR variations at the disk wall have time scales comparable to 
that of the ejection of HH objects. However, the disk instabilities near the location of the 
dust sublimation zone cannot induce accretion events on time scales shorter than the length 
of time required to accrete matter from the wall to the star, and this is also measured in 
thousands of years or longer. 

While it seems unlikely that changes in the region of the DSZ could directly impact the 
accretional activity on time scales shorter than centuries, changes in the inner gaseous disk 
will have much smaller time scales and could impact both the DSZ and the accretion rate 
(an subsequent gas outflow in the jets). In the case of HD 163296, a “rapid” (i.e., time scales 
of weeks to months) change in the scale height of the inner gas disk could lead directly to 
changes in the illumination of the dust wall. This will create a widening zone of subliming 
grains whose outer radius will move outward from the star until the grains are cool enough 
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to survive. The changes observed in the dust scale height (whose rapidity will be on the 
order of t z ) and i nn er radius (governed by how rapidly the grains can sublime) are consistent 
with such a scenario. 


4-3.2. Magneto- Rotational Instabilities 

If the inner dusty disk is partially ionized, then a magneto-rotational instability (Balbus 
k Hawley 1991a, b) could also develop. Such instabilities have been recently invoked to 
describe gas accretion without equally efficient dust accretion in actively-accreting objects 
with “cleared” inner disk holes (Chiang k Murray-Clay 2007; Kretke k Lin 2007). Such 
an application is particularly relevant to disks where the main inner disk wall is far beyond 
the location of the DSZ for refractory grains, but may coincide with that of volatile grains 
(the “snow line”). Here the instability is invoked to help explain the presence of accretion 
in what would otherwise be non-accreting systems. The main requirements are ionizing 
radiation from the star - usually assumed to be X-rays from both the accretion columns 
and stellar corona with magnetic fields embedded in the ionized gas. Both HD 31648 and 
HD 163296 possess hot accreting gas (and HD 163296 is an X-ray source), and the highly 
collimated outflows suggest significant magnetic fields in the star-disk region, providing the 
main ingredients for MRI processes. In the models presented here, the X-rays and ultraviolet 
emission from the accretion shock are responsible for some of the heating of the disk. Changes 
in accretion onto the star must alter the heating of the exposed dust wall. While the MRI 
model has generally been invoked to provide a mechanism for steady accretion, the details 
of how time-dependent changes would occur, and their possible effect on the DSZ, are not 
known. But such a mechanism might provide a direct link between activity near the surface 
of the star (perhaps related to jet activity and HH object production) and changes in the 
DSZ. 


4.3.3. The X- Wind 

Both accretional heating from within the disk and high energy stellar and accretion 
shock photons photons hitting the surface can cause the gaseous disk to become ionized. In 
the X-Wind models of Shu et al. (1994) (and see the comments in Shu et al. (2007)) the 
stellar magnetic field is connected to that of the differentially-rotation gas disk. The model 
not only provides a mechanism for the production and collimation of the outflowing jets 
observed in the stars discussed here, but opens the possibility that reconnection events in 
the magnetic field could alter both the outflow rate and structure in the ionized gas disk at 
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the DSZ interface. Should the X-ray and ultraviolet flux from the accretion shock ionize the 
surface of the DSZ, then it may also be affected. 


4- 3. 4- Planetary Perturbations 

Finally, we are observing these stars during a phase of disk development when planets 
could be forming. Planet formation has been suggested for the truncation of the thick dusty 
disks far beyond the DSZ for refractory grains in the disk of TW Hya (Calvet et al. 2002), 
and the thinning of material interior to it, effectively creating a class of “transition” disk 
system beyond the Meeus Group I and II classification. If the structure of the inner wall is 
in any way affected by such a body having an orbit close to the current sublimation distance 
(about 0.24 AU), then its orbital period would be on the order of months. So far, we have 
no indication of changes on time scales that small in HD 163296, and this time scale is an 
order of magnitude shorter than that of the ejection of the HH objects. 


4.4. Other Observable Consequences 

If changes are occurring in the location and height of the DSZ, these should be observable 
in a number of different ways. If some of the stellar radiation intercepts the outer flared disk, 
then changes in scale height of the inner rim will alter the degree of shadowing of the outer 
disk. This will produce changes in both the scattered light and the far-IR emission. Multi- 
epoch observations of the scattered light coronagraphic images of HD 163296 suggest such 
changes might have been detected (Wisnewski et al. 2007). Far-IR monitoring with sufficient 
time coverage and precision does not exist for this star, but such variations may have been 
seen in the thermal emission of the HAeBe star SV Cep (Juhasz et al. 2007). 

If the DSZ undergoes changes in exposure to stellar radiation, there would be two other 
important consequences that could be tested. As exposure increases, newly exposed dust 
will find itself above the sublimation temperature, and the DSZ must retreat from the star 
(Figure 16). When and if the gas opacity and/or scale height is re-established, some of the 
vaporized material, if not driven away by stellar or disk winds, may re-condense, and the 
DSZ would gradually fill back in. This might be observable with NIR interferometers, but 
would require regular monitoring programs. 

Because the emissivity of grains is selectively poorer for smaller grains than larger ones, 
the smaller grains would be the first to vaporize. Any surviving grains would be large and 
would undergo annealing if their temperatures were greater than about 800 K (Hallenbeck 
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et al. 1998, 2000; Hill et al. 2001). Once irradiation of the disk returns to lower values, 
grain condensation, occurring over time spans of months, could add additional crystalline 
material to this zone. Ablation of the dust in this region by stellar winds and photon 
pressure would gradually redistribute some of this material further out, but would still be 
preferentially be concentrated at its source, the DSZ region. Because the ratio of grain 
mass to surface are increases with increasing grain size (at least for compact grains), these 
processes will transport the smaller grains more effectively, leaving the larger grains closer to 
the DSZ. Infrared spectroscopy of the silicate band emission using the Very Large Telescope 
Interferometer (van Boekel et al. 2004) indicates that the inner few AU of these stars are in 
fact characterized by both larger grains and more crystalline material, relative to the outer 
regions, as would be expected with this scenario. Multi-epoch spectroscopic monitoring of a 
few of these sources with the VLTI would be useful investigating this possible mechanism. 


4.5. One Size Does Not Fit All: The Remarkable Case of DG Tau 

While migration of the DSZ may be sufficient to explain the variations seen in the SEDs 
of HD 31648 and HD 163296, this mechanism cannot explain the behavior of other PMS 
objects exhibiting near-IR variability. 

Wooden et al. (2000) first reported unusual changes in the 8-13 pm spectrum of the 
T Tauri star DG Tau. In 1996, its IR spectrum consisted of a very weak silicate emission 
band, with spectral structure consistent with the presence of both amorphous and crystalline 
silicates. A year later, the emission band was weaker (almost nonexistent), and the amor- 
phous component was undetectable. In early 2001, a dramatic change was captured by two 
different instruments: DG Tau exhibited an an overall increase in 3-5 pm emission, and the 
onset of a silicate band in absorption (Woodward et al. 2004). The 2001 BASS spectrum, 
shown in Figure 17, is (so far) unique for a PMS star in showing the presence of crystalline 
silicate material in absorption. Between 2002 and 2004, the star returned to its previous 
“conventional” behavior. But in 2006 it underwent another 3-5 pm outburst, but on this 
occasion, the silicate band was in emission! The breadth of the emission band requires mean 
grain sizes in excess of a few pm, smaller than those present in its quiescent state, but larger 
than the ones responsible for the 2001 absorption event. 

The relatively simple DSZ expansion and migration scenario is clearly inapplicable to 
DG Tau. Rather, a scenario more like the “ablation” one discussed by Vinkovic & Jurkic 
(2007) seems more likely. This picture is consistent with the observation by Kitamure et al. 
(1996) of the expansion of the disk in DG Tau, which they surmise may be driven by a stellar 
wind. However, while the expansion derived from their CO radial velocity measurements 
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are primarily confined to the outer disk, the lofting observed in 2006 requires warm grains, 
indicating that this radial outflow of disk material may originate near the DSZ in this star. 


5. Conclusions 

We have presented observations of two “isolated” Herbig Ae stars that exhibit significant 
variability in their infrared emission. Our immediate finding is that any model of the SED 
based on data obtained in a single epoch cannot describe the nature of source completely. 
In this sense, there is no “the model” for any such object, and this problem is exacerbated 
further when modeling includes the use of disparate data sets obtained on different epochs 
(often using diverse photometric systems). Our goal is to provide the sort of multiwavelength 
multiepoch data that models require if they are going to be capable of explaining all of the 
data in a self-consistent manner with the fewest changes in model parameters. 

HD 31648 and HD 163296 test our ability to understand these disks, because both 
exhibit features that are characteristic of PMS disks systems that are both very young and 
simultaneously relatively old in terms of evolutionary development. Both are sources of polar 
jets of outflowing material and HH objects which must be supplied by the active accretion 
of gaseous material close to or onto the star. That is, both also exhibit other indicators of 
highly active inner regions which are associated with accretion shocks, such as far-UV and 
X-ray emission and high-temperature ionic gas species. Yet their SEDs in the mid to far 
infrared have led them to be classified as stars where the outer disk has undergone significant 
grain coagulation and settling to the point where much of this region is shadowed by the 
innermost disk, which has been forced to expand vertically by a paucity of accreting gas. 

In both stars, interferometric observations (PTI and IOTA) indicate the presence of 
an inner disk similar to that expected in models that include a puffed-up inner wall due to 
direct illumination by the central star. Geometries like these are consistent with the Group 
II classification and the faintness of the coronagraphic observations of the scattered light in 
the outer disks. Nevertheless, the sum of the data present us with a number of challenges. 

For HD 31648, Isella et al. (2006) have remarked that the predicted SED derived from 
models of the interferometry measurements is “only roughly consistent with the photometric 
values.” Their models generally predict an increase in brightness (in units of AF>) with 
wavelength that is not a good match to the photometrically-determined SED. In their model 
SEDs, the predicted brightness at 3.5 pm is consistently too large. An examination of the 
dates when the PTI observations were made (Eisner et al. 2004) indicate that they were 
obtained over a time span of a year, during which time HD 31648 may have been varying in 
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brightness at the wavelength that these observations were made. The nearly “flat” 1-2 /xm 
spectral shape seen in the models derived from the IOTA observations (Monnier et al. 2006) 
are consistent with the photometric data only if HD 31648 was in a bright state when it was 
observed. To properly derive the structure of these disks, nearly simultaneous photometric 
and interferometric measurements will be required, and the effects of hot accreting gas 
included. 

HD 163296 presents a different problem. Here, the interferometry data, the disk shadow- 
ing, and the presence of a strong 1-5 gm hump are all consistent with the vertically expanded 
inner wall at the DSZ. But how can the disk be illuminated directly by the star as required 
by such passive disk models, yet still be actively accreting gas at a rate large enough to drive 
the bipolar gas flows and HH objects? This requires that the equatorial optical depth of the 
inner gas disk be less than unity, or its vertical height be much smaller than the height of 
the dust sublimation wall. This apparent dichotomy may be solved if the inner gas disk is 
much thinner than the disk wall, and the accretion rate small, but sufficient to shadow at 
least part of the dust disk mid-plane, as in the model of Muzerolle et al. (2004). We may be 
observing many of these sources just prior to them reaching the transition disk phase, where 
the accretion of gas is diminishing, but not yet terminated, and may be sporadic. Should 
the scale height of the disk have a vertical “step” - going from the inner region goverened 
by gas opacity to the outer one governed by dust opacity - it might result in an inner dust 
sublimation zone of complex structure, where the equatorial regions are indeed shadowed by 
the gas, but the upper regions exposed to unattenuated stellar radiation. 

In both stars we may be witnessing the final phase of active accretion. HD 31648 has 
one of the lowest detected accretion rates in intermediate-mass PMS stars (Grady et al. 
2008) and is also amongst the oldest of these stars to exhibit active accretion. The time 
scale over which changes in the innermost gas disk and jet activity could be on the order of 
a few years, while that of the DSZ itself must probable be longer. If so, repeated exposure 
and shadowing of the DSZ will lead to a region inside the DSZ where grains are repeatedly 
heated and cooled. It is here that one would expect to be able to anneal amorphous grain 
material into a crystalline phase if they survive sublimation, or condense them directly into 
crystalline form if they do not. Repeated heating and cooling dust material in the DSZ 
would enhance the production of crystalline material there, as has been observed by VLTI 
spectroscopy. 

In the models by Isella et al. (2006) for the inner disk rim, grain sizes “larger than ~ 
1.2 /im and possibly much larger” are required to explain the interferometry and the gross 
features of the SEDs of the Herbig stars they studied. However, much larger grains would 
be inconsistent with the grain sizes needed to produce a strong 10-pm silicate band in the 



27- 


“outburst” of HD 163296, if that material is associated with the region of the inner disk 
wall. 

Part of the difficulty with matching the models derived from interferometry measure- 
ments with the observed SEDs may result from intrinsic variability of the emission. In- 
terferometric measurements made at different epochs may yield different results, and such 
changes have been recently confirmed by Tannirkulam k Monnier (private communication). 
Complete u-v plane coverage needs to be obtained quickly for a single epoch, along with 
simultaneous SED photometry, before such observations can be effectively used to construct 
realistic models of the inner disk region. On the positive side, multi-epoch coverage will 
provide a very powerful tool for studying physics of these regions. 

Finally, the case of DG Tau serves as a warning that the actual structure of the inner 
regions of these systems must be complex and highly variable. Expecting a simple time- 
independent model to describe this system is doomed to failure. Models that include even 
the most rudimentary dynamical changes are needed. 

We are far from being able to claim that we truly understand the structure of PMS 
disks. For a better understanding of these systems, it is crucial that repeated observations 
with NIR interferometers be obtained on a more regular basis over the course of at least a few 
years, accompanied by nearly-simultaneous photometry. It is our hope that measurements 
such as these will provide useful constraints on future modeling of the inner disk regions of 
these objects. 
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A. Appendix A: Infrared Photometry & Flux Calibration 

While filter technology and our understanding of atmospheric transmission has gradually 
led to the evolution of “standard” filter sets for near infrared (wavelengths shorter than 5 
/zm) measurements (Simons & Tokunaga 2002; Milone & Young 2005), the situation in the 
8-14 /mi telluric transmission window is less standardized. 

A variety of infrared magnitude systems have also been in used during the past 3 decades. 
In many cases, these derive from the different filter systems being used, so that transforma- 
tions between systems is required. In general, these may be relatively precisely calibrated 
[for example those between six different systems described in detail by Glass (1985)]. Some 
are based on a system where a Lyr has a magnitude of 0.0 at all wavelengths, while in others 
it is the color index of an “average” A0 star that is defined as being 0.0, while a Lyr itself 
may be anywhere between +0.02 and -0.02 mag (Glass 1985) or even brighter. 

Somewhat more puzzling are systematic differences between some filter systems and 
higher spectral resolution spectrophotometric systems. For example, a CMa is a favorite 
calibration star because it is bright, and because its infrared spectrum contains relatively 
weak photospheric lines, making the calibrations insensitive to the exact wavelengths of 
the detectors. Nevertheless, many photometric systems use magnitudes for this star in 
the 2-14 /zm region that are roughly 1.50 mag brighter than a Lyr (Gehrz et al. 1974), 
while spectrophotometric measurements made with BASS (Russell & Mazuk 1998) and other 
systems (Cohen et al. 1995) indicate a value closer to 1.35 mag., a 15% difference. 

When actually comparing real absolute flux values (instead of magnitudes), as is usual 
for SEDs, the choice of an absolute flux calibration system is important. These generally fall 
into three broad classes: those based on the absolute flux of the Sun, those based on stellar 
atmosphere models, and those based on terrestrial blackbody sources. The precise system 
chosen will lead to differences in the derived fluxes. For example, the absolute flux of a Lyr 
at 8.7 /mi on the BASS system in units of A F\ is 1.90xl0 -11 W m 2 , while that of Gehrz et 
al. (2005) is 1.8xl0 -11 W m 2 , and that of Cohen et al. (1992) is 1.71xl0 -11 W m 2 , a range of 
10%. A comparison of the BASS and Cohen systems, along with that used by Sitko (1981) 
for the first multiwavelength SED study of the stars in this paper, are shown in Figure 16. 

New work has recently been shed light on the application of stellar photospheric models 
to a Lyr. It has been known for some time that a Lyr sits above the main sequence, and 
that the width of its spectral lines implied a small value of the projected rotational velocity 
v sin i. Optical interferometry by Peterson et al. (2006) support the notion that the star is 
actually flattened due to rapid rotation (93% breakup speed) and seen nearly pole-on. They 
show that their model predicts excess emission over the standard spherical stellar atmosphere 
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models by an amount that grows from 1.6% at 0.67 pm to 7.2% at 3.69 pm. 

Because of all of the effects described above, the blind use of different relative and 
absolute flux systems can introduce spurious flux differences of up to 17%. 

Despite these differences, for the purpose of studying the variability of SEDs, the precise 
system to be used is not as important as using (as much as possible) the same system for 
all observations. For the purpose of this study, we have used the measurements of Russell 
& Mazuk (1998) to define the relative brightness of our calibration stars in all systems, and 
have used the BASS absolute flux calibration throughout. We have re-examined the data 
reductions used for the observations of the HAeBes reported previously by Sitko (1981) and 
Sitko et al. (1994), with the aid of the original observing logs, in order to make a better 
estimation of the true errors and eliminate possible systematic effects due to changes in 
seeing, sky transparency, and calibrations stars used, as indicated in those logs. The dates 
of the observations, telescopes, and instruments used for the infrared observations used here 
are listed in Table 1. 

In addition to the infrared observations listed in Table 1, we have also used the Low 
Resolution SpectroPhotometer (LRSP) observations obtained with the 0.4-m reflector of the 
Pine Bluff Observatory of the University of Wisconsin. These data are the same ones used 
by Sitko (1981) for determining the first mutliwavelength SEDs of these two stars. The 
instrument utilized a movable grating to sample the spectrum between 3390 and 5840 A, 
using a photomultiplier tube and photon-counting system. The spectral resolution was set 
by an aperture mask that could sample the spectrum in bandpasses of 40 and 400 A. For the 
observations presented here, the spectrum was sampled at 17 wavelengths using the 40 A 
bandpass mask. The main use of these observations is to tie the infrared data to UV obser- 
vations obtained with the International Ultraviolet Explorer. They have a distinct advantage 
over filter photometry since they are capable of providing a “clean” Balmer discontinuity not 
possible with the standard broad-band U filter. These data also usually include the peak flux 
of the stars, which are important for subsequent modeling calculations. They also allow a 
comparison with filter photometry from the literature, in order to detect possible variability 
at visible wavelengths. 
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HD 163296 



A (^m) 

Fig. 1. — The spectral energy distribution of HD 163296. From shortest wavelengths to 
longest: solid lines - IUE merged spectrum (LWR05626 & SWP06566 on 1979 Sep 19) and 
Pine Bluff LRSP (1979 Jun 22, 1979 Aug 12, & 1979 Sep 17); open squares between 0.36 /an 
and 4 /xm - UBVRI photometry [15 sets over 9 nights, from de Winter et al. (2001)], JHKLM 
photometry [4 nights between 1980 & 1986, from de Winter et al. (2001)], and 2MASS; solid 
line between 13 and 40 /tm - Spitzer IRS (2004 Aug 28); filled circles - BASS (1996 Oct 14); 
open squares from 12 to 100 /xm - IRAS photometry; open squares longward of 100 /xm - 
millimeter & submillimeter photometry from Henning et al. (1994), Henning et al. (1998), 
Mannings & Emerson (1994), and Isella et al. (2007). 
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Fig. 2.— The spectral energy distribution of HD 31648. From shortest wavelengths to 
longest: solid lines - IUE merged spectrum (LWR04278 & SWP04940 on 1979 Apr 15) 
and Pine Bluff LRSP (1979 Mar 10); open squares between 0.36 /im and 4 pm - UBVRI 
photometry [1998 Oct 25 from Oudmaijer et al. (2001)], JHK photometry [1998 Oct 28, 
from Eiroa et al. (2001)], and 2MASS; filled circles - BASS (1996 Oct 14); open squares 
from 12 to 100 pm - IRAS photometry; open squares longward of 100 pm - millimeter &. 
submillimeter photometry from Mannings et al. (1997), Mannings & Sargent (1997), and 
Pietu et al. (2006). 
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HD 163296 
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Fig. 3. — The infrared flux of HD 163296 on 14 dates, organized as 10 separate epochs. Due 
to the proximity in time and similarity in flux level, the following can be considered single 
epochs: 1996 Oct 10 (ISO) & 14 (IRTF), 1998 May 12 (IRTF) & 16 (CST), 2002 Jul 18 
(Lick) & 29 (IRTF), 2005 Jul 06 (Lick) & 26 (IRTF). The ESO data are from de Winter 
et al. (2001); the CST data are from Eiroa et al. (2001). The KPNO data were previously 
published in Sitko (1981). 
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HD 163296 



Fig. 4. — The combined VNIRIS+BASS spectroscopy of HD 163296 for 2002 and 2005. The 
lowest curve and diamond symbols are the difference between the flux levels measured by 
VNIRIS and BASS, respectively, for the two epochs. The thicker smooth solid line is a model 
spectrum of a mixture of graphite and silicate grains (the latter being Mgo.8oFeo.2oSi0 4 ). The 
particle size distribution used was that of Dust Impact Detection System (DIDSY) in situ 
measurements of Comet lP/Halley (McDonnell et al. 1987, 1991), with size cutoffs at 0.03 
and 0.25 /zm. Also shown in the figure is the difference in the 3-13 /im flux of HD 31648 
between 1996 and 2004 (triangles). The thinner smooth solid line is the same grain model 
as for HD 163296, but with a temperature of 1600 K, a minimum grain size of 0.25 /zm 
[consistent with the model of Isella et al. (2006)] and a maximum size of 3.0 /im. Such a 
large minimum size cannot explain the changes in flux for HD 31648. Unless the grain size 
range is unrealistically restricted, a “large” grain population cannot account for the relatively 
strong presence of the silicate band in the differenced spectra of these two systems. 
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HD 163296 
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Fig. 5.— Light curve of HD 163296 at 3.77 pm. The diamonds are ESO photometry, while 
the filled circles are derived from BASS/IRTF spectrophotometry. The UT dates of the 
observations are given as yymmdd. 
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Fig. 6. — The infrared flux of HD 31648 on 10 dates. The 1979 KPNO data were previously 
published in Sitko (1981). The KPNO and MLOF data between 1980 and 1984 are previously 
unpublished. 
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HD 31648 & HD 163296 



Fig. 7 — Variability of HD 31648 and HD 163296 at 3.77 pm versus versus 10.5 pm. With 
the exception of the 2002 outburst, the flux near 10 pm is relatively constant in HD 163296, 
even when the 3 pm emission changes. For HD 31648, changes at 3 pm have so far been 
accompanied by similar changes at 10 pm. 
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HD 1632968 Dustfiles 
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Fig. 8.— The spatial distribution of the grain files used for the model of HD 163296. The 
envelope includes small “ISM” grains, while the disk contains the “medium” sized grains in 
its outer layers, and “large” grains closer to the disk mid-plane. The transition from the 
blue to yellow indicates where the density reaches the “mid-plane” regime (number density 
of H 2 > 10 10 cm -3 ). The model for HD 31648 was similar, except that the model fit was 
somewhat better without the transition to the largest grains in the mid-palne. 
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Fig. 9. — The temperature and density structure for the model of HD 163296, at three 
different radial scales The color indicates how the temperature and density progress from 
lowest to highest values, with blue the lowest, followed by red, orange, yellow, and finally 
white for the highest. The “truncated envelope” has a density approximately ten orders of 
magnitude lower than the disk midplane, but receives direct illumination from the star. 
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Fig. 10. — A model of the IR emission of HD 31648 using the Monte Carlo radiative transfer 
code of Whitney et al. (2003a, b, 2004), and described in the text . For this model the BASS 
spectrum from 1996 has been used. For simplicity, we show the model results for only three 
disk inclinations, 0° (upper solid line), 90° (lower solid line), and 30° (intermediate solid 
line). The last of these is close to the inclination values of Simon et al. (2001) and Pietu et 
al. (2006). For the stellar photospheric emission, a Kurucz model atmosphere with T e ff = 
8250 K and log g = 4.5 was used. 
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Fig. 11 — Models for HD 31648 on two epochs. The 3-13 fj,m data in the triangles consists 
of the BASS spectrum from 1996, with shorter-wavelength 1998 observations of similar flux 
levels (lower open squares), from the EXPORT project (Eiroa et al. 2001). The dotted 
line is the model for the 1996-98 data that is shown in Figure 10. The shorter wavelength 
observations represented by diamonds were taken from the 2MASS catalog, and are also not 
contemporaneous with the BASS observations, but of approximately the same flux state. 
The filled circles are the BASS data from 2004. For the 2004 epoch, the fit was achieved 
with the same model as for the earlier epoch, except the disk scale height was increased from 
0.0076 AU to 0.0130 AU, and had a very slight increase in envelope mass. The upper squares 
are the JHK observations from Oct. 1980 at MLOF (see Appendix A and Table 1), are not 
contemporaneous, but used to illustrate the general shape of the SED at a higher flux state. 
The dashed line is the same photospheric model used in Figure 10. These particular models 
slightly underpredict the strength of the silicate band, and suggests that the very surface 
layers of the disk may have grain sizes slightly smaller than those used here. 
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HD 163296 2005 
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Fig. 12. — A model of the IR emission of HD 163296 in 2005. The upper and lower solid lines 
are for inclinations of 0° and 90°, respectively. The middle curve is 50°, consistent, within 
the uncertainties, of the values reported in the literature (Mannings & Sargent 1997; Grady 
et al. 2000; Wassell et al. 2006; Isella et al. 2007). For the stellar photospheric emission, a 
Kurucz model atmosphere with T e // = 8750 K and log g = 3.5 was used. 
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Fig. 13. — A model of the IR emission of HD 163296 in 2002 using the same model as 2005, 
but increasing the disk scale height from 0.012 to 0.018 AU. Also shown is the same model, 
with the added “outburst” component discussed in the text, and shown in Figures 4 and 14. 
The BASS and VNIRIS data are the same as those shown in Figure 4. 
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HD 163296 2002 
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Fig. 14. — A model of the IR emission of HD 163296 in 2002, using the additional component 
responsible for the outburst has been added to the 2005 fit. 
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Fig. 15. — The Ca II triplet, Pa j3, and Br 7 lines on the two epochs for which VNIRIS 
data are available. The equivalent width of the Ca II feature at 0.850 nm was approximately 
4.7 A in 2002 and 6.6 A in 2005. The decreased strength in 2002, when the hump was 
strongest, suggests that the change may have been due to either increased occultation of the 
triplet-emitting region as the hump outburst proceeded, or by the addition of an unresolved 
absorption component in the line profiles. No significant changes were observed in the 
strengths of the hydrogen lines. 
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Fig. 16. — Pictogram of possible changes occurring in the inner disk region of HD 163296, 
as described in the text. The shapes of these regions are meant to be illustrative of one 
way that the structure might change, not an exact representation of the structure, since the 
current radiative transfer code does not include shadowing by the inner gas disk. 
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Fig. 17 — DG Tau, observed on 4 dates between 2001 and 2006. Also shown is one spectrum 
obtained using the SL1 grating of the Infrared Spectrograph (IRS) of the Spitzer Space 
Telescope (from the SST Archive). 
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Fig. 18. — A comparison of the various flux calibration systems. The data show the absolute 
flux of a Lyr for: BASS (upper curve), Cohen et al. (1992) (lower curve), earlier Sitko 
(1981) (open squares), and empirical values from Campins, et al. (1985) and Rieke et al. 
(1985) (filled circles). 
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Table 1. Observations 


Star UT Date Telescope Instrument 


HD 31648 


HD 163296 


31 Aug 1979 
03 Jan 1980 
06 Oct 1980 
08 Oct 1980 
02 Dec 1984 
14 Oct 1996 

05 Aug 2004 

11 Dec 2006 

1 Sep 1979 
02 Sep 1979 
14 Oct 1996 

12 May 1998 

18 July 2002 
29 July 2002 
29 May 2004 

06 July 2005 
26 July 2005 

19 May 2006 


KPNO 1.3m 
KPNO 1.3m 
MLOF 1.5m 
MLOF 1.5m 
KPNO 1.3m 
IRTF 
IRTF 
IRTF 

KPNO 1.3m 
KPNO 1.3m 
IRTF 
IRTF 
Lick 3m 
IRTF 

MLOF 1.5m 
Lick 3m 
IRTF 
IRTF 


Bolo 

Bolo 

Uplooker3 

Uplooker3 

Hermann 

BASS 

BASS 

BASS 

Bolo 

Otto 

BASS 

BASS 

NIRIS 

BASS 

BASS 

VNIRIS 

BASS 

BASS 


Note. — “Bolo” refers to the bolometer of the Kitt Peak 
National Observatory (KPNO). The system utilized one “fla- 
vor” of JHKLMNQ broad-band filters, as well as intermediate- 
bandpass filters located in the 10 pm region, useful for deter- 
mining the strength of the silicate emission band. “Otto” 
and “Hermann” at KPNO were InSb photometers with sim- 
ilar JHKLM filters. “Uplooker3” was yet another InSb sys- 
tem used at the Mount Lemmon Observing Facility (MLOF), 
which at that time was operated by the University of Cal- 
ifornia at San Diego and the University of Minnesota, and 
is currently supported by The Aerospace Corporation and 
the University of Minnesota. For all of these observations, 
detailed records of the sky conditions have been preserved, 
and in many cases (i.e. the KPNO bolometer data) raw data 
output suitable for re-processing have survived the quarter- 
century of time since the data were obtained. 



Table 2. Model Parameters 


Parameter HD 31648 HD 163296 


Stax T e ff (K) 

Distance (pc) 

A v (mag) 

Disk Midplane Dust File 3 * * 
Disk Dust File 3 
Envelope Dust File 3 
Disk Outer Radius (AU) 

Disk Inner Radius (AU) 

Disk Scale Height - Inner (AU) b 
A - Disk Mass Density Exponent 0 
B - Disk Scale Height Exponent** 

Disk Accretion Rate (Mq yr -1 ) 

Disk Mass (Mg) 

Envelope Max Radius (AU) 

Envelope Min Radius (AU) 

Inner Cavity Wall Opening Angle (deg) 
Envelope Mass Infall Rate (Mg. yr _1 ) c 
Envelope Mass (Mg yr _1 ) c 
Accretion Shock Luminosity Fraction 
C - Surface Brightness Exponent* 3 


8250 

8750 

122 

140 

0.05 

0.10 

Medium 

Large 

Medium 

Small 

Small 

Small 

250 

450 

0.196/0.245 

0.29/0.35 

0.0076/0.0130 

0.012/0.018 

1.6 

1.99 

0.6 

0.99 

4.0xl0- 9 

8.0x10-® 

0.8 

0.4 

10 

10 

0.24 

0.83 

10 

20 

0xl0~ 7 /1.2xl0' 6 

7-OxlO -7 

,0xl0- a /6.4xl0- 8 

2.6x10'® 

1.1% 

11.0% 

-3.1 

-3.1 


Note. — The sublimation temperature for all the models was set to be 
1600 K. According to Gail (2002), for disk mid-pane conditions, the sublima- 
tion/condensation temperatures for enstatite (MgSiOs) and forsterite (Mg 2 Si 04 ) 
silicates are between 1300 and 1500 K. For some other minerals that could provide 
hot dust emission, like perovskite (CaTiOa), this temperature can be significantly 
higher. 

3 The “small”, “medium”, and “large” grains refer to the interstellar medium 

grain model of Kim et al. (1994), the HH30 scattered light grain model of Cotera 
et al. (2001), and the HH30 disk thermal emission model of Wood et al. (2002), 
respectively. 

b A range in values are for the “low” and “high" flux states. 

“Describes the radial mass density distribution within the disk.: p oc R~ A 

d Describes the change in scale height of the disk with radial distance: H oc R B . 

B = 1 corresponds to a disk with a conical opening. B > 1 indicates a concave 
upward flared disk, while B < 1 is for a convex upward “anti-flared” disk. 

e The predicted A=1.14pm surface brightness distribution, where the surface 
brightness <xR _c . The value listed for HD 163296 matches the observed value, 
based on coronagraphic imaging with the Hubble Space Telescope. Although not 
detected in the earlier survey of Grady et al. (2005), the disk of HD 31648 has 
recently been detected (Stecklum et al. 2007), and is steeper than predicted by 
the current model. 






